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Abstract 



In this paper we investigate the effect of the resummation of threshold logs 
on the rate for B — > X s ^. We calculate the differential rate dT/dE^ including 
the infinite set of terms of the form a™log n+1 (l — x) and a™log n (l — x) in 
the Sudakov exponent. The resummation is potentially important since these 
logs turn into log(2E cu t /mi,) , when the rate is integrated from the lower cut 
x = 2E cu t/rrib to 1. The resummed rate is then convolved with models for the 
structure function to study whether or not the logs will be enhanced due to 
the fermi motion of the heavy quark. A detailed discussion of the accuracy of 
the calculation with and without the inclusion of the non-perturbative effects 
dictated by the B meson structure function is given. We also investigate the 
first moment with respect to (1 — x), which can be used to measure A and Ai. 
It is shown that there are some two loop corrections which are just as large 
as the a 2 s (3o term, which are usually expected to dominate. We conclude that, 
for the present energy cut, the threshold logs do not form a dominant sub- 
series and therefore their resummation is unnecessary. Thus, the prospects 
for predicting the rate for B — ► X s ^ accurately, given the present energy cut, 
are promising. 
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I. INTRODUCTION 



The process B — > X s "y is considered fertile ground for discovering physics beyond the 
standard model due to the fact that its leading order contribution is a one loop effect. Of 
course, any hope of finding new physics is predicated on our ability to control the theoretical 
prediction within the standard model. This decay rate can be calculated in a systematic 
expansion in a s , l/(M w ,m t ) and l/mj,. Tremendous effort has gone into calculating the 
decay rate at next to leading order in the strong coupling. The calculation is broken into 
three stages. First the full standard model is matched onto the four fermion theory. It is 
then run down to the scale m&, after which the inclusive decay rate may be calculated in 
heavy quark effective field theory IJ. 

It is in this last step that certain effects arise which jeopardize the accuracy of the 
calculation. In particular, due to the fact that there is an imposed energy cut on the 
photon, a new small parameter enters the calculation, namely, the distance to threshold 
1 — 2.£7 7 /m&. As i±7 7 approaches m&/2, we begin to probe larger and larger distances, until 
finally, uncontrollable non-perturbative effects dominate. 

The subject of this paper will center on the threshold logs, log(l — 2E cut /rrib), which 
become parametrically large as E cut approaches the end point. As is well known, these logs 
appear because the limited phase space obstructs the KLN cancellation of the IR sensitivity. 
Thus for large enough values of E cut the perturbative expansion breaks down, and must be 
resummed. The leading double logs resum into the well known Sudakov exponent. However, 
the resummation beyond leading log accuracy is much more complicated. The next to 
leading logs were resummed in Ref. || in a systematic fashionF] utilizing the moment space 

1 We take the opportunity to correct several typos in Q. After Eq. (40) the definition of B should 
include a factor of Cp. In Eq. (57) both 7 and B should have plus signs in front of them and in 
Eq. (74) gfW should be replaced by 7W. In Eq. (75) the argument of the exponential should be 
log Ngi +g 2 . 
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factorization ||. However, since the result was given in moment space, an analysis of 
the result in terms of measurable quantities is difficult. Indeed, in Ref. [0] it was found 
that in moment space the sub-leading logs in the exponent actually dominate the leading 
logs. Thus, it is important that we make sure that the physical prediction not suffer the 
same consequences. Previously in the literature the subleading logs in x space (x being the 
variable 2E 1 /rrib) were naively estimated by exponentiating the one loop subleading log HH. 
In Ref. [H an attempt was made at including the subleading logs in a non-trivial fashion, 
however, it was not shown that all the logs of a given order were summed. 

Here, we will build upon the results in Ref. || and present the resummed differential 
rate at next to leading log accuracy in x space. We then integrate the rate as a function 
of the energy cut to get the rate to be compared with experiment. Then the effects of 
the fermi motion, which in principle could enhance the threshold logs, are included in this 
resummation. Finally, the net theoretical errors incurred due the the unknown shape of the 
structure function responsible for the fermi motion are estimated. 



II. THE ONE LOOP RESULT 

By combining the operator product expansion (OPE) and the heavy quark effective 
theory (HQET), it is possible to calculate the decay spectrum for inclusive heavy meson 
decays in a systematic expansion in a s and AqcdA 7 ^ 0- I n the endpoint region, E~ — > 
rrib/2, both the perturbative and AqcdA^ expansions break down. The non-perturbative 
corrections can be resummed into a structure function H,0], which will be discussed in 



Section VI 



The calculation of inclusive B decay rates begins with the low-energy effective Hamilto- 
nian 

AC 8 

H ef f = --^^Ec^o.Gu), (I) 

j= i 

where Gp is the Fermi constant, are elements of the Cabibbo-Kobayashi-Maskawa matrix, 



3 



Ci(fj) are Wilson coefficients evaluated at a subtraction point /i, and Oj(/i) are dimension 
six operators. For B — ► X s 7, the only operators that give a relevant contribution are 

g 

7 = j^mbSLaCr^F^bRa, (2) 

Here e is the electromagnetic coupling, g is the strong coupling, mj is the 6 quark mass, 
F^ u is the electromagnetic field strength tensor, G% u is the strong interaction field strength 
tensor, and T a is a color SU(3) generator. Near the endpoint, the decay rate is dominated 
by the O7 operator. In particular, the decay rate due to the O7 operator is 



l_dT 

r dx 

where x = 2E 1 /mi > and 



(2x 2 - 3x- 6)x + 2(x 2 - 3) log(l - x) 
^ ~ 3(1 - x) 

X<1 V / 



(3) 



Gl\V*V tb \ 2 aC*ml 
10 = 32^ • (4) 

As x — > 1 this contribution diverges as log(l — x)/(l — x). The integrated rate is finite due 
to virtual corrections at x = 1. Only decays with large photon energies can be detected, 
due to large background cuts with the current experimental cut on the photon energy given 



by £L > 2.1 GeV [10 . 



III. THE SYSTEMATICS OF THE EXPANSION 



The perturbative series near the endpoint schematically takes the form 

dT logfl — x) 1 

— = C n a s — hCi 2 a s - hCi 3 a s 

dx 1 — x 1 — x 

.,, 2 log 2 (l-x) 2 log(l-x) 
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All non-integrable functions here are tacitly assumed to be "plus" distributions. When we 
integrate form 5 to 1, then each power of log n_1 (l — — x) turns into log n (l — 5). Given 
this series we may ask, how close to the endpoint can we study the spectrum and still expect 
to get the right answer? Clearly, if a s log 2 (l - 5) C 1 we may simply truncate the series at 
the first term. To go closer to the endpoint, i.e. larger values of 5, it would seem that we 
must sum the complete lower triangle along a diagonal, a daunting task. For instance, if we 
just sum the first column, then we can not let a s log 2 (l — 5) ~ 1 since at some point the 
terms in the next column may grow. However, this naive criteria is incorrect. The reason 
for this is that the series is known to exponentiate into a particular form. In fact, the series 
resums into a function of the form 

log[r(5)] = log(l - 8) 9l [a s log(l - 5)} + g 2 [a s log(l - 5)} + a s g 3 [a s log(l - 5)] + ■ ■ ■ . (6) 

This result implies that the aforementioned triangle assumes a definite structure. The re- 
summation does not sum the entire triangle. There will always be cross terms between 
higher order and lower order terms in the exponent, which arise in the expansion of the 
exponential, that have not been kept. What the resummed form does tell you though, is 
that these cross terms can be neglected. The important point to notice about this structure 
is that if we truncate the series in Eq. (D, and assume that the last term we kept is 0(1), 
then there is no "growth" in higher orders in a s as there is in the case of organizing the 
calculation in terms of the columns. 

If we keep the first term in the expansion of g\ we get the usual Sudakov double logarithm. 
In that case S must satisfy the condition a 2 log 3 (l — 5) < 1. Note, that this does not allow 
log 2 (l-<5) to become arbitrarily large (practically this will not be an issue). In general, 
then if we expand gi up to 0(a™ _1 ), the requirement becomes [a s log(l — 5)} n log(l — 5) < 1. 
Thus, once a s log(l — 5) approaches one, we must for all intensive purposes include the entire 
gi, as well as g 2 . 

On a practical note it is clear that we can not let the logs become arbitrarily large. 
As S approaches one we begin to probe momenta of order Aqcd, where the perturbative 

5 



approach stops making sense. Equivalently, we are asking questions about details of the 
hadronization process once we reach the resonance regime. Presently, the experimental 
limit on 5 is 2(2.1)/m& — 0.88.Q Thus a s log 2 (l — 5) ~ .95, if we do not include the factor of 
1/tt which accompanies each power of a s . If we include the factor of it, then it may well be 
that we have not resummed the dominant piece of the series, and resumming non-logarithmic 
terms could be just as important. The question of the inclusion of a factor of 7r is a dicey 
numerical issue given the size of coefficients in the expansion. Thus it is necessary to perform 
the resummation to determine whether or not the logs form a dominant sub-series. 

IV. THE RESUMMATION 

The resummation is performed in moment space, where the rate factorizes into a short 
distance hard part (H^), a soft part (SV) and a jet function (Jjv). Using renormalization 



group techniques it is possible to show that [11 



a N J N = exp [\ogNgi(x) + flfcOe)] , (7) 
where x — a s( m ft)A)l°giV, and 

9i = -7T^-[(1 " 2X) log(l - 2 X ) - 2(1 - X ) log(l - X)}, (8) 
2npoX 



log(l-2x)-21og(l-x) 



— — - log(l - x) ~ TT-r log(l - 2x) 

4vr^ 27t/? 



+ Ilog 2 (l-2 X )-log 2 (l- X ) 
+ ^ [log(l - 2 X ) - log(l - x)] • (9) 

In the above, f3 = (11 CU - 2N f )/(12ir), fa = (17 -5C A N f - 3C F N f )/(2Art 2 ), and 
= C A (67/18 - tt 2 /6) - 10 Tr JV>/9. 



2 Here x is defined in terms of the quark mass, this will be modified when we include the effects 
of the structure function. 
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To go back to x-space, we must take the inverse-Mellin transform of Eq. (|7|) which is 
given by 

HV 1 rC+ioo 

T x = 2^i c MNX ' NdN (10) 

Then we use the identity JT^ ] 

— / x -iV e io6(Ar)F[ as iog(iV)] diV= _ " ( Q ( 1 _ x ) e ^K0 x \ l + E ( a Al) 
27ri Jc-ioo ax L J 



where 

E(a s ,/) = E«"E^ (12) 
fc=i j=o 

represents subleading log contributions and 

/ = -log(-log(x)). (13) 

It is important at this point to note a crucial difference between this calculation and re- 
summations carried out in threshold production. In threshold production (Drell-Yan for 
instance) the overall sign of g\ is flipped (in standard schemes). This arises as a consequence 
of the fact that the Sudakov suppression in the parton distribution function overwhelms the 
analogous suppression in the hard scattering amplitude. Thus, after subtraction, the overall 
sign of the exponent changes. This difference in sign changes the nature of the inverse Mellin 
transform. In particular, note that in the case of a negative exponent (our case) the inverse 
Mellin transform is integrable even if we ignore the derivative acting on the step function, 
which is multiplying a function that is manifestly zero. Whereas in the case of the positive 
sign, the 0-function is crucial to define the inverse Mellin transform in a distribution sense. 
The expansion of the series for the Drell-Yan case has a non-integrable pole on the positive 
axis in the Borel plane. Which is to say that using the above approximation has introduce 
subleading terms which nonetheless become numerically important due to spurious factorial 
growth. Thus, as pointed out by Catani et. al. to avoid any large corrections/ambiguities 
one should not use Eq. ([TT|) but instead one must perform the inverse Mellin transform 



exactly JT2|. In our case, the expansion of the analytic result also leads to a factorially 
divergent series, but it is sign alternating and therefore Borel summable. More simply put, 
the x space form is integrable. 

To calculate the inverse Mellin transform with next-to-leading log accuracy {i.e. sum all 
logs of the form a n s log n (l — 5) in the exponent) we will perform the integral numerically. In 
so doing we we choose the integration contour such that the constant C above is chosen to 
lie to the left of the Landau pole singularity and to the right of all other potential poles. 



V. EFFECTS ON THE TOTAL RATE AND FIRST MOMENT 

Let us now consider the effects of the resummation on the partonic rate. We will start 
by expanding the function g\ in a s and systematically improving the approximation by 
including more terms. The expansion of g-y leads to a convergent series with a radius of 
expansion given by x max = 1 — exp[— l/(2a s /?o)]. To integrate the rate from 5, we use the 
fact that the first moment is one, and thus instead integrate from to 5, thus avoiding the 
region where our approximation of the moments breaks down. If we keep up to 0(a™) in 
this expansion, then the terms we drop are 0(a" +1 log n+2 (l — 5)). Thus keeping more terms 
in the expansion allows us to take a s log(l — 5) closer to one. Figure 1 shows the percentage 
contribution to the total rate stemming from the resummation, as a function of the photonic 
energy cut. We see that as we include more terms in the expansion of g\ in the exponent, 
the expansion converges to a fixed rate. Also, in Fig. 1, we show the rate which includes 
all the subleading logs of the form [alog(l — 5)] n in the exponent, which is calculated by 
keeping the full form of g\ and g 2 and performing the inverse Mellin transform numerically. 
We see that overall the effect of resummation is small. Moreover, the subleading corrections 
are of the same size as the leading corrections. This is NOT because the expansion is ill 
behaved, but because the subleading terms are not truly suppressed compare to the terms of 
order a™ log n+1 (l — 5), simply because the logs is not that large. We do not have a dominant 
sub-series to sum. The inclusion of g 2 gives a larger result because the 0(a s ) piece in g 2 has 
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FIG. 1. The effects of resummation on the rate, integrated from 5 to 1. The solid curve is the 
leading log result. The sparse dotted curve is the numerically inverted next to leading log result. 
The tight dotted curve is the rate expanding g\ to 0{a s ). The long dashed curve is expanding g± 
to 0(a 2 s ), the dot-dashed curve is expanding g\ to O(af), and the short dashed curve is expanding 



a larger coefficient than the 0(a s ) piece in g\. 

Thus, the resummation is effectively just including some piece of the two loop result. 
This can be clearly seen in Fig. 2 where we show the resummed double log rate with the 
0(a s ) piece subtracted out, and the 0{a 2 s ) piece from the expansion of the double log 
resummation. We see that the results are nearly identical. Also in Fig. 2 we show the 0(a 2 ) 
derived from expanding the full result, including g 1 and g 2 , performing the inverse Mellin 
transform exactly using 



We see that the 0{a 2 s ) piece coming from formally subleading log terms is actually larger 



gi to 0(a A s ). 



| rC+ioo 
2%i JC-ioo 



x~ N log" (N)dN 




(14) 
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FIG. 2. The solid curve is the resummed double log rate, with the 0(a s ) piece subtracted out. 
The dotted curve is the 0(a 2 ) piece of the resummed double log rate. The dashed curve is the 0{a 2 s ) 
piece of the full resummed rate, including both g\ and gi . 

than those coming from the leading logs terms. This was hinted at in the results of the 
moment space calculation in Ref. ||. However, this is not because the expansion for the 
rate is converging poorly, as the effect is still small compared to the lower order contribution, 
but because the logs are just not that large. 

We may use also our results to pick out certain terms in the two loop calculation. In 
particular, we may expand the exponent and determine the coefficients of the terms of 
order a 2 log 3 (l — x)/{l — x) and a^log 2 (l — x)/{l — x). Lower order logs will not be 
correctly reproduced given the fact that we have dropped terms of order a s in the exponent. 
Expanding Eq. (^) to 0(a 2 ) we find the terms 



dT 

dx 



8a 2 \og 2 (l-x) 



0(a2) 



+ _M 2/%7r + — 



14 \ logfl - x) 



(15) 



9tT 2 1 — X 7T 2 V 37 1 — X 

Recently, the BLM correction to the first moment of dT/dx has been calculated. This 
correction, proportional to a 2 f3o, typically is a good approximation to the full two loop 
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result. The interest in calculating the corrections to the first moment is to get a better 
handle on the extraction of the HQET parameters A and Ai. If we compare the piece 
proportional to /3q in Eq. ([15]), we find agreement with the corresponding piece found in 
Ref. [|13j0. However, we also note that the contribution of the second term in Eq. fllS]), not 
proportional to (3$, is just as big as the BLM corrections calculated in Ref. ]T3| . This does 
not rule out the possibility of a cancellation with other terms such that the BLM result in 
the two loop result still dominates. However, it does put into question the issue of BLM 
dominance, and calls out for the complete two loop calculation of the first moment. 



VI. CONVOLUTION WITH THE STRUCTURE FUNCTION 

As is well known, as one probes the spectrum closer the end point, the OPE breaks 
down, and the leading twist non-perturbative corrections must be resummed into the B 
meson structure function Formally we may write the light cone distribution function 

for the heavy quark inside the meson as 

f(k+) = (B(v) | h v S(k+ - iD + )h v | B(v)). (16) 

While the shape of this function is unknown, the first few moments of f(k + ), 

A n = J dk + k n + f(k + ) 

= (B(v) | hJ(iD + ) n h v | B(v)), (17) 

are known; A — 1, A\ — 0, and A 2 = — Ai/3. f(k+) has support over the range — oo < 
k + < A. It has been asserted that the support below —A dies off exponentially 
formal proof has been given. 

The effects of the fermi motion of the heavy quark can be included by convoluting the 
above structure function with the differential rate, 



, but no 



5 Our definition differs for (3q differs by a factor of 4ir from [13] 
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/A dk + f(k + )^(mt), (18) 



-m b dEj 

where dYp/dE^ is the rate Eq. (||), or the resummed version Eq. fllCf ), written as a function 
of the "effective mass" ml = m b + k + , i.e., x = 2E 1 /ml = MEXs/fnl- The new differential 
rate is now a function of xb = 2E 1 /Mb- The addition of the structure function resums the 
leading- twist corrections and moves the endpoint of the spectrum from E 1 = m b /2 to the 
physical endpoint E 1 = The cut rate may then be written as 



r 



H 



IE, 



cut 



M 



B 



A 

2E cu t—'^b 



IE. 



cut 



(19) 



where T p [2E cut / (m b + k + )} is the partonic rate with a cut at x p = 2E cut /{m b + k + ). From 
this result we can see that at smaller values of k + the infra-red logs will be enhanced because 
the phase space available to gluon emission is curtailed. On the other hand, the structure 
function will be suppressed at values of k + < —A. Thus, we expect any enhancement to be 
vitiated by the effects of the structure function. 

We will use following two different ansatze for the shape of the structure function ||: 

Fi{k+) = iVi(l - A; + /A) a e (1+a)fc+/A , (20) 
F 2 (k+) = N 2 (l - k + /A) b e- c{1 ~ k+/A)2 . (21) 

The parameters N{, a, b, and c can be determined from the known moments v4o,i,2- Note that 
both of these forms have exponential suppression for k + < —A, which need not be true. The 
values of A and Ai have been extracted from the data and are highly correlated. We will use 
the central values (A = 0.39, Ai = —0.19) and the one sigma values (A = 0.28, Ai = —0.09) 
and (A = 0.50, A x = -0.09) determined in flj . 



In Fig. 3 we show the rate as a function of the lower cut on the photon energy, for the 
ansatz in Eq. ( p0|) and with the central values for (A, Ai). In this plot we also show the rate 
with and without perturbative resummation without the inclusion of the structure function, 



The true end point of course takes into account the final state masses. 
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FIG. 3. The rate integrated from E cu t- The solid curve is the resummed rate convoluted with 
the structure function in Eq. $2(\), using the central values of (A, Ai). The shaded region shows the 
uncertainty due to varying the values of (A, Ai) by one sigma. The dotted curve is the one loop rate 
without the structure function. The long dashed curve is the resummed rate without the structure 
function. The dot- dashed curve is the one loop rate convoluted with the structure function, using 
the central values of (A, Ai). The short dashed curve is the one loop rate and the 0{a 2 s ) piece from 
the resummation, convoluted with the structure function. 

and the one loop rate as well as the one loop rate plus the 0(a 2 ) piece from the expansion 
of the resummed rate, convoluted with the structure function. As expected, the effects of 
both perturbative and non-perturbative resummation are enhanced as the photonic energy 
cut is increased. We see again that the dominant piece of the resummation is coming from 
the 0{a 2 s ) term. Furthermore, the inclusion of the structure function does not significantly 
enhance the effects of the resummation for these choices of structure functions. If they died 
off more slowly for large negative k + , we would expect the logs to be enhanced. The shaded 
region is the resummed rate using the one sigma values for (A, Ai), with again the structure 
function defined in Eq. (|2"0|) . We see that, as A is increased, the effects of resummation 
become more important, as one would expect since the width of the primordial distribution 
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FIG. 4. Rate convoluted with the different structure functions. The solid curve is the rate 



convoluted with the Eq. (21), and the dashed curve is the rate convoluted with Eq. (21) 



is increasing. Varying the width (A, Ai), is one way of determining the uncertainty due to 
our ignorance of the structure function. Another method would be varying the functional 
form of the structure function itself {i.e. changing the higher moments). 

In Fig. 4, we compare the rates for the two differing functional forms given by Eq. (|20|) 
and Eq. (pip. We see that the rate is not as sensitive to the higher moments as it is to 
the lower moments (A, Ai). Note that for these particular choices of the structure function 
the fourth and fifth moments differ by only O(50%). So perhaps a more thorough search 
of the possible form of the structure functions should be explored. Assuming that the 
structure function is well behaved (read physically motivated), the dominant uncertainty 
in the rate due to our ignorance of the structure function can be removed once we have a 
better determination of (A, Ai). 

Finally, we may investigate the convergence of the expansion of g\ in the exponent with 
the inclusion of the structure function. In Fig. 5 we show the rate for the series expansion 
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0.65 




FIG. 5. Rate convoluted with the structure function Eg. The dotted curve corresponds to 

expanding g\ to 0(a s ), the dashed curve is for g\ expanded to 0(a 2 s ), and the dot-dashed curve is 
expanded to 0(oP s ). The solid curve is the rate with the full g\. 

of gi up to 0(al), we see that it is indeed well behaved and converges rapidly. This again 
is signaling that the resummation is not collecting a set of dominant terms. 

VII. CONCLUSIONS 

We conclude that resummation is not necessary, and does not increase the accuracy 
of the prediction, when the photonic cut is near 2.1 GeV. Therefore, for the purpose of 
calculating the total rate, it is consistent to convolve the structure function with the 0(a s ) 
partonic rate. The dominant errors will then come from our ignorance of the structure 
function and higher order, in a s , corrections to the rate, as well as uncertainties in T due 
to the dependence on the b quark mass (i.e. A). We estimate that the uncertainties due 
to the structure function are at the 20% level, though we believe that a more thorough 
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search of the space of structure functions should be performed. If we assume that the size 
of the 0(af) corrections are typically of the size of the pieces of the two loop result which 
we pick off from our resummed rate, then the uncertainty coming from these corrections is 
on the 10% level. On the other hand, if we wish to extract the structure function from the 
measurement of the spectrum for b — > sj, in order to utilize it in the extraction of V u b, the 
resummation will most probably be necessary. 
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